Charmonium production in heavy-ion collisions from the SPS to LHC 
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Using the two-component model that includes charmonium production from both initial nucleon- 
nucleon hard scattering and regeneration in the produced quark-giuon plasma, we study J/ip pro- 
duction in heavy-ion collisions at the SPS, RHIC and LHC. For the expansion dynamics of produced 
hot dense matter, we use a schematic viscous hydrodynamic model with the specific shear viscosity 
in the quark-gluon plasma and the hadronic matter taken, respectively, to be twice and ten times 
the lower bound of 1/An suggested by the AdS/CFT correspondence. For the initial dissociation and 
the subsequent thermal decay of charmonia in the hot dense matter, we use the screened Cornell po- 
tential to describe the properties of charmonia and perturbative QCD to calculate their dissociation 
cross sections. Including regeneration of charmonia in the quark-gluon plasma via a kinetic equation 
with in-medium chamonium decay widths, we obtain a good description of measured J/i/; nuclear 
modification factors in Pb-|-Pb collisions at ^/snn ~ 1.73 GeV at SPS and in Au-f An collisions at 
y^SNN = 200 GeV at RHIC. A reasonable description of the measured nuclear modification factor 
of high transverse momenta J/^ in Pb-f Pb collisions at y^sjvjv = 2.76 TeV at LHC is also obtained. 

PACS numbers: 



I. INTRODUCTION 



Since J/tp suppression was first suggested by Matsui 
and Satz as a signature of quark-gluon plasma (QGP) for- 
mation in relativistic heavy- ion collisions [l|, there have 
been many experimental [J, Q and theoretical studies 0- 
|8|l on this very interesting phenomenon; see, e.g., Refs. 
[3, [13] for a recent review. The original idea of Matsui 
and Satz was that the color screening in the produced 
QGP would prohibit the binding of charm and anticharm 
quarks into the J/ip and thus suppress its production. 
However, lattice QCD calculations of the J/^ spectral 
function have since shown that the J/tJj can survive 
above the critical temperature for the QGP phase tran- 
sition [nl, As a result, the study of J/tp suppression 
in relativistic heavy-ion collisions has been changed from 
being a signature of the QGP to a probe of its proper- 
ties. Indeed, we have recently shown in a two-component 
model, which includes J/ip production from both initial 
hard nuclcon-nuclcon scattering and regeneration from 
charm and anticharm quarks in the produced QGP, that 
the in-medium effect on J/tJj interactions in the QGP 
can affect the J/ip nuclear modification factor and ellip- 
tic flow in Au-|-Au collisions at y^SNN = 200 GeV at 
the Relativistic Heavy Ion Colhdcr (RHIC) [Hi. In the 
present study, we extend this study to J/ip production in 
Pb-I-Pb collisions at the higher energy o f W snn = 2.76 
TeV at the Large Hadron Collider (LHC) [ll[l| and also 
at the lower energy of y^SNN — 17.3 GeV at the Super 
Proton Synchrotron (SPS) [1]. Furthermore, a schematic 



'Electronic address: lsongtsoo@yonsei.ac.kr 
t Electronic address: khan@co mp.tamu.edu| 
■f Electronic address: |ko@comp.tamu.edu| 



viscous hydrodynamic model is used to include the effect 
of viscosity on the expansion dynamics of the produced 
hot dense matter that was neglected in our previous stud- 
ies. We find that the two-component model can give a 
good description of the experimental data from heavy- ion 
collisions at these different energies. 

To make the present paper self contained, we briefly 
review in Sec. [H] the two-component model for J/tp pro- 
duction, in Sec.[ni|the schematic causal viscous hydrody- 
namical model used in modeling the expansion dynam- 
ics of produced hot dense matter, and in Sec. IIVI the 
in-medium dissociation temperatures and thermal decay 
widths of charmonia. Results obtained from our study for 
the J/^J nuclear modification factors in heavy-ion colli- 
sions at SPS, RHIC and LHC are then presented in Sec. 
fVl Finally, a summary is given in Sec. IVII 



II. THE TWO-COMPONENT MODEL 

The two-component model for J/ip production in 
heavy-ion collisions [l6l [l7j includes contributions from 
both initial hard nucleon-nucleon scattering and regener- 
ation from charm and anticharm quarks in the produced 
QGP. For initially produced J/tp's, their number is pro- 
portional to the number of binary collisions between nu- 
clcons in the two colliding nuclei. Whether these J/tp's 
can survive after the collisions depends on many effects 
from both the initial cold nuclear matter and the final 
hot partonic and hadronic matters. The cold nuclear 
matter effects include the Cronin effect of gluon-nucleon 
scattering before the production of the primordial J/ip 
from the giuon-giuon fusion [l^; the shadowing effect 
due to the modification of the gluon distribution in a 
heavy nucleus |l9|: and the nuclear absorption by the 
passing nucleons |20l422| . In our previous work [23| on 
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J /ip production in heavy- ion collisions at RHIC, we have 
considered only the most important nuclear absorption 
effect. In this case, the survival probability of a primor- 
dial J /ijj after the nuclear absorption is given by [l^, [25| 



'S'cnm(b, S 



TAB{h s) / ^^^^'/'A(s,2;)ps(b- s,z') 
xexpj - - 1) j dz^p^(s,ZA)CTabs| 
xexpj - (B -I) dzBPB{h - s,zs)crabs|, (1) 



where b is the impact parameter and s is the transverse 
vector from the center of nucleus A; TabO^, s) is the nu- 
clear overlap function; pAi^, z) is the density distribution 
in the nucleus; (Tabs is the J/ip absorption cross section 
by a nucleon. For the latter, it is obtained from p+A col- 
lisions and has values of 4.18 and 2.8 mb for the SPS and 
RHIC, respectively 0, H^. Presently, there arc no p-|~A 
data available from the LHC. Since the cross section for 
J/ijj absorption is expected to decrease with increasing 
energy [2l| , we consider in the present study the two ex- 
treme values of and 2.8 mb to study its effect on the 
J/ip yield in heavy-ion collisions at LHC. 

Although the shadowing effect has usually been ne- 
glected in heavy-ion collisions at SPS and RHIC, this 
may not be justified at LHC. In the present study, we 
thus include also the shadowing effect for heavy-ion col- 
hsions at LHC using the EPS09 package [l3 . The shad- 
owing effect is expressed by the ratio Rf of the parton 
distribution ff^{x,Q) in a nucleus to that in a nucleon 
jnucioonj-^ Qj niultiplicd by the mass number A of the 
nucleus, i.e.. 



Rfix,Q) 



AH 



clcon 



(2) 



In the above, x = niT / \/snn, with niT being the trans- 
verse energy of the produced charmonium and ^Jsnn 
being the center-of-mass energy of colliding nucleons, is 
the momentum fraction and Q = rriT is the momentum 
scale. Assuming the shadowing effect is proportional to 
the path length, we can then express the spatial depen- 
dence of Rf^ as S m in : 



i?f(s,x,Q) 



1 ^ j^J dzpA{s,z) 



Rf{x,Q)-l JdzpA{0,z)' 



(3) 



where is a normalization factor determined from the 
condition 



dzpA{s, z)Rf{s, X, Q) = Rf{x, Q). (4) 



The shadowing effect reduces the survival probability 
of a primordial J /ip after the nuclear absorption (Eq.(IT])) 
by the factor Rf{s, x, Q)Rg {h — s, x, Q). Taking the mo- 
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4.2 MeV to be the average J ftp trans- 
1.96 TeV fH, we obtain the 
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0.897 


TO (fm/c) 


1.0 


0.9 [43] 


1.05 [43] 




ri/s 


0.16 


0.16 [43] 


0.2 [43] 





TABLE I: Parameters for J jip production and the firecylinder 
expansion in Pb-fPb collisions at ^Jsnn = 17.3 GeV at SPS 
and at ^snn = 2.76 TeV at LHC and in Au-|-Au collisions 
at y^SNN = 17.3 GeV at RHIC. da^^^^Jdy and da^^/dy are, 
respectively, the differential J/^ and cc production cross sec- 
tions in rapidity in p-l-p collisions; fxa, /i/.'(2S), and are, 
respectively, the fraction of J/ip production from the decay 
of Xc, i^' , and bottom liadrons in p-|-p collisions; R^ is the 
gluon shadowing effect on charm production; and tq and r]/s 
are the thermalization time and the specific viscosity of the 
produced QGP. Also shown in the last column are the pa- 
rameters for the feed-down contribution to the production of 
J/tp's of transverse momentum pr > 6.5 GeV at LHC. 



value of the ratio R^^{x,Q) given in Table |T] for charm 
production in heavy-ion collisions at the LHC. 

For the hot partonic and hadronic matter effect, the 
model includes the dissociation of charmonia in the QGP 
of temperatures higher than the dissociation temperature 
and the thermal decay of survived charmonia through 
interactions with thermal partons in the expanding hot 
dense mater. Since the number of produced charm quarks 
in relativistic heavy-ion collisions is not small, charmo- 
nia can also be regenerated from charm and anticharm 
quarks in the QGP. The effect of thermal dissociation and 
regeneration of charmonia on the number Ni of charmo- 
nium of type i is taken into account via the rate equation 

m 



dNj 
dr 



-T,{N.,-NP), 



(5) 



where r is the longitudinal proper time, while N^'^ and 
Ti are, respectively, the equilibrium number and ther- 
mal decay width of charmonia and will be discussed in 
SeclTVl 

Since charm quarks are not expected to be completely 
thermalized either chemically or kinetically during the 
expansion of the hot dense matter, the fugacity param- 
eter 7 and the relaxation factor R are introduced to de- 
scribe their distributions. Assuming that the number of 
charm and anticharm quark pairs does not change dur- 
ing the fireball expansion, the fugacity is obtained from 
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where is the number of cc pairs produced in an 

A+B colhsion; tIq and Uh are, respectively, the number 
densities of open- and hidden-charm hadrons in grand 
canonical ensemble; V is the volume of the hot dense 
matter; and Iq and Ii are modified Bessel functions re- 
sulting from the canonical suppression of charm quarks 
in heavy- ion collisions [UH^- For the relaxation factor, 
it is defined as i?(r) = 1 — exp[— (t — ro)/req] with the 
relaxation time Ton = 3 fm/c of charm quarks in the QGP 
taken from Ref. [7| and tq being the initial thermaliza- 
tion time. 

Since charmonia can only be regenerated in the QGP 
of temperature below the dissociation temperature Ti, 
the number of equilibrated charmonium of type i in the 
QGP is 



j^R fVe{T, - T), 



(7) 



where is its number density in grandcanonical ensem- 
ble; / is the fraction of QGP in the mixed phase and is 
1 in the QGP; and 6{Ti — T) is the step function. 

For the initial charmonium number Ni and the charm 
quark pair number Ncc^ they are obtained from multiply- 
ing their respective differential cross sections in rapidity 
daf^/dy and da^^/dy in p-l-p collisions [l^,!!!,!!^! by the 
number of binary collisions A^coU in heavy-ion collisions. 
Since only the J/i/j production cross section at y/s — 7 
TeV dH has been measured in p-l-p collisions at LHC, 
its value at ^/sWn — 2.76 TeV is obtained by using a 
linear function in ^/s to interpolate from the measured 
values at ^/s = 1.96 TeV by the CDF Collaboration at 
the Fermi Lab ^ to the one at = 7 TeV at LHC. 
The cross section for cc pair production at ^/snn = 2.76 
TeV is then determined by assuming that the ratio be- 
tween the J/tJj and cc pair production cross sections is 
the same as that at RHIC. In Table U we list the differ- 
ential cross sections for J/ip and cc pair production in 
p-|-p collision at SPS, RHIC and LHC that are used in 
the present study. 

Since J/ip production in p+p collisions includes the 
contribution from the decay of excited charmonium 
states, the cross section da^f^^Jdy shown in Table U is 

the sum of the production cross sections for the J / ip and 
its excited states. For p-l~p collisions at SPS, we use the 
global average values of the fractions /^^ = 25% from the 
Xc decay and /^'(2S) = 8% from the ^p' decay [s^. The 
cross sections for J/ip, Xc and ip' production in a p-l-p 
collision at the SPS are then given, respectively, by 
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where o-j^^ is the cross section for J/ tp production with- 
out the feed-down contribution, and 'Br' denotes the 
branching ratio. 

For p-|-p colisions at RHIC, the fractions of J/ip^s from 
Xc and ip'{2S) decays are taken to be f^^ = 32% and 
ftpi = 9.6 %, respectively, based on recent experimen- 
tal results by the PHENIX Collaboration ^ . Since the 
fractions of J/t/i's from Xc and ^'{2S) decays are not 
known at LHC, we use the values inferred from p+p anni- 
hilation at y/s^N = 1.96 TeV by the CDF Collaboration 
at the Fermi Lab. It was found in these reactions that 
among promptly produced J/ip's, about 64% are directly 
produced and about 29.7% from the Xc decay, and both 
are approximately independent of the J/ip transverse mo- 
mentum [ssj . This leads to the fraction of promptly pro- 
duced J/ip's from the ip'{2S) decay to be 6.3%. Using 
the experimental result that promptly produced J/ijj^s 
constitute about 89% of measured J/'0's [2^, we ob- 
tain the fractions /^^ = 26.4% and /^'(2S) = 5.6% of 
measured J/tp^s that are from Xc and ip'{2S) decays, 
respectively. Since the fraction of prompt J/i/j's is re- 
duced to 79% for J/ip of transverse momentum pT > 6.5 
GeV [1^, the fractions of measured J/t/i's of Pt > 6.5 
GeV that are from Xc and ip'[2S) decays are reduced to 
fy^^ = 23.5% and /^'(2S) = 5.0%, respectively. Besides 
the contribution from excited charmonia, the decay of 
bottom hadrons can also contribute to J /ip production 
in high-energy collisions. This contribution increases sig- 
nificaiitly with px as shown in the experiments by the 
CDF [il, CMS m, LHCb ^ and ATLAS Col- 
laborations. The fraction is between 5 % and 10 % for 
Pt < 3 GeV, depending on the rapidity of the J /ip^ then 
increases to more than 40 % at pt ~ 15 GeV, and reaches 
60-70 % for Pt above 25 GeV. On the average, about 11 % 
of produced J ftp's are from the decay of bottom hadrons 
in p -I- p annihilation at ^/s = 1.96 TeV at the Fermi Lab 
[29| . and the fraction increases to about 21 % for J ftp's of 
transverse momentum px > 6.5 GeV [2^. These values 
and other input parameters used in the present study are 
shown in Table ID 



III. A SCHEMATIC VISCOUS 
HYDRODYNAMIC MODEL 

For the expansion dynamics of the hot dense matter 
formed in relativistic heavy-ion collisions, we describe it 
by a schematic causal viscous hydrodynamic model re- 
cently developed in Ref. . It is based on the assump- 
tion that all thermal quantities such as the energy den- 
sity, temperature, entropy density, and pressure as well as 
the azimuthal and space-time rapidity components of the 
shear tensor are uniform along the transverse direction 
in the hot dense matter. Assuming the boost-invariance 
and using the (r, r, 0, r/) coordinate system 



BvitP' -> JftP + X)' 



(8) 



1 , t + z 
— In , 

2 t- z' 
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r = v^x^ + y'^i (/i = tan ^(y/x), 



(9) 



then the following equations are obtained from the usual 
Israel-Stewart viscous hydrodynamic equations: 



(10) 



ar(^(7.))-^4(<+<?), (11) 



A 



(12) 



(13) 



In the above, T"^"^ = (e + Prju"^ — Pr is the time- 
component of the energy-momentum tensor, tt^ = t^tt'^'^ 

and ttJj' = t^tt'''' are, respectively, the azimuthal and the 
space-time rapidity component of the shear tensor; rjs 
and T,r are the shear viscosity of the hot dense matter and 
the relaxation time for the particle distributions, respec- 
tively; e = ^drirjr) + ^dr{rVr-fr) with = l/\/l " 

in terms of the radial velocity Vr', A = ttR^ with R being 
the transverse radius of the uniform matter; and (• • •) 
denotes average over the transverse area. For the ra- 
dial flow velocity that is a linear function of the radial 
distance from the center, i.e., jr^r = lRR{'f/R), where 
R = dR/dr and = - R^, we have (7^) = 

l+7|i?V2, il^vf.) = 7i^V2, (7.) = 2(7i-l)/(37ii?^), 
and {jrVr/f) — jrR/R. With the energy density e and 
pressure p related by the equation of state of the matter 
through its temperature T, Eas. (fT0|) - (|T3| are four simul- 



taneous equations for T, R, and ttJ^', and can be solved 
numerically by rewriting them as difference equations. 

For the equation of state of the produced dense matter, 
we use the quasiparticle model with three flavors for the 
QGP phase [4l|, \^ and the resonance gas model for the 
HG phase. As to the specific shear viscosity rjs/s, where 
s is the entropy density, its value in the QGP is taken to 
be 0.16 for SPS and RHIC, and 0.2 for LHC [H, while it 
has the same value of 5/2tt in the HG The specific 
viscosity in the mixed phase is assumed to be their linear 
combination, i.e., (v/ s)qgp f+{v/ s)hg{^- f), where / is 
the fraction of QGP in the mixed phase. The initial ther- 
malization time is taken to be 1.0 fm/c for SPS, which 
has usually been used, and 0.9 fm/c and 1.05 fm/c for 
RHIC and LHC, respectively Although the initial 

thcrmalization time for RHIC is 0.6 fm/c in ideal hydro- 
dynamics [45I , the nonzero viscosity generates additional 



transverse flow |23| and requires a late thcrmalization to 
fit the experimental data on pT spectra and elliptic flows. 
This is the same reason for the later thcrmalization at 
LHC in viscous hydrodynamics. 

The initial local temperature of produced matter can 
be calculated from the equation of state and the local 
entropy density, which we parameterize as [4^ |46| 



dr] 



(1 - a)^-^ + a ncoii 



(14) 



with a = 0, 0.11 and 0.15 for SPS, RHIC and LHC, 
respectively [46l - l48| . The number density fipart(coii) in 
Eq. dH]) is defined as AiVpaj.t(coii)/(T'o^a^^2/), where 
^A'part(coii) is the number of participants (binary col- 
lisions) in the volume tqAxAij of the transverse area 
AxAy and is obtained from the Glauber model with the 
inelastic nucleon-nucleon cross sections of 30, 42 and 64 
mb for SPS, RHIC and LHC, respectively The 
factor C is determined by fitting the multiplicity of final 
charged particles after the hydrodynamical evolution to 
the measured one. 

Assuming the same chemical freeze out temperature 
Tf = 160 McV for all charged particles, their pseudora- 
pidity distribution at midrapidity is then [Tsj 



dN, 



ch 



dr] 



■^"^Di J dpT Pt J rdr 



rriTi^ ""dydpT 
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Xin 



Pt sinh p 




rriTi cosh p 


[ Tf \ 




[ Tf \ 



, (15) 



where p = ta,nh.~^ {vr) ■ The summation i includes all 
mesons lighter than 1.5 GeV and all baryons lighter than 
2.0 GeV. In including the contribution from the decays 
of particles, we simply multiply their pseudorapidity dis- 
tributions by the product Di of their decay branching 
ratio and the number of charged particles resulting from 
the decay. We have thus neglected the difference between 
the rapidity of the daughter particles and that of the de- 
cay particle. Also, we have used the thermal momentum 
distributions at chemical freeze out as well as during the 
expansion of the hot dense matter, thus ignoring the vis- 
cous effect on the particle momentum distributions as it is 
only important for particles of large momenta jsif . From 
the multiplicities of charged particles per half participant, 
((iA^ch/d?7)/(A^part/2), which are roughly 2, 4 and 8.4 in 
central collisions of Pb4-Pb at y^s^N = 17.3 GcV at SPS, 
of Au+Au collisions at ^/Tnn — 200 GeV at RHIC, and 
of Pb-HPb at y^sWN = 2.76 TcV at LHC, respectively 
[4^[53, we obtain the corresponding values of 14.6, 18.7 
and 27.0 for the parameter C in Eq. ([H)) . 

In Fig. [1] we show the temperature profile along radial 
direction at initial thcrmalization time in heavy-ion col- 
lisions at SPS, RHIC, and LHC from the viscous hydro- 
dynamics. Defining the firccylindcr as the region where 
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FIG. 1: Temperature profiles along the radial direction at 
initial thermalization time as functions of radial distance in 
central Pb+Pb collisions at y^sjvjv = 17.3 GeV at SPS and 
^/sNM ~ 2.76 TeV at LHC, and in central Au+Au collisions 
at ^/sNN = 200 GeV at RHIC from viscous hydrodynamics. 



FIG. 2: Average temperatures of firecylinder as functions of 
time in central Pb+Pb collisions at y^sjviv = 17.3 GeV at 
SPS and ^/s^ = 2.76 TeV at LHC, and in central Au+Au 
collisions at ^sjvjv = 200 GeV at RHIC from the viscous 
hydrodynamics . 



the initial temperature is above Tc = 170 MeV, its trans- 
verse radius in the case of viscous hydrodynamics has 
values of 6.5, 6.6 and 7.1 fm in central collisions at SPS, 
RHIC and LHC, respectively. The time evolution of the 
average temperature of the firecylinder determined from 
the schematic hydrodynamic model is shown in Fig. [2] 
The initial average temperatures at SPS and RHIC are 
218 and 269 MeV, respectively, and are consistent with 
those extracted from the experimental data on dileptons 
at SPS m and on direct photons at RHIC [H. The pre- 
dicted initial average temperature in heavy-ion collisions 
at LHC is 311 MeV. 

For non-central heavy-ion collisions where the ini- 
tial geometry of the transverse area is an ellipse, the 
schematic viscous hydrodynamic model described here 
needs to be extended. For simplicity, the present model 
is used by taking the circular transverse area to be the 
same as that of the ellipse as in Ref. [42| based on a 
parameterized firecylinder model. 



IV. THERMAL PROPERTIES OF CHARMONIA 

To describe the properties of charmonia in QGP, we 
need the potential between heavy quark and its anti- 
quark at finite temperature. Although some information 
on this can be obtained from the lattice gauge theory 
(55I [56| or from the static QCD [s^, HI], we use in the 
present study the extended Cornell model that includes 
the Debye screening effect on color charges [H^. The 
Cornell model [g^l was devised to imitate the asymptotic 
freedom and confinement of the QCD interaction with 
a Coulomb-like potential for short distance and a linear 
potential for long distance. In the QGP, the linear poten- 



tial becomes weaker due to the Debye screening between 
color charges, leading to the screened Cornell potential 

in 



V{r,T) 



-M(T)r 



(16) 



with a = 0.192 GeV^ and a = 0.471. The screening mass 
/i(r) depends on temperature and is given in thermal 
pQCD by 



(17) 



where Nc is the number of colors, Nf is the number of 
light quark flavors, and g is the QCD coupling constant. 
In the limit of /i — >■ 0, we recover the original Cornell 
potential. 

The wavefunctions and binding energies of charmonia 
in the QGP are obtained by solving the Schrodinger equa- 
tion with the screened Cornell potential. With the bind- 
ing energy Eq defined as [5^ 



eo = 2mc H E, 



(18) 



where the charm quark mass is taken to be rric = 
1.32 GeV and E is the eigenvalue of the Schrodinger 
equation, we show in Fig. |3] the binding energy of J/^p 
as a function of temperature for the case of g — 1.87. 
It is seen that the J/^p becomes unbound or dissociated 
in the QGP for temperatures above ~ 300 MeV. As in- 
dicated by Eq. ([T7| . the J/ip dissociation temperature 
decreases as the QCD coupling constant g increases. This 
is shown in Fig. 0] not only for the J/ip but also for its 
excited states Xc and ip' . In obtaining the dissociation 
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FIG. 3: Binding energy of J /il) in the QGP as a function of 
temperature for the QCD couphng constant g = 1.87. 



temperatures for Xc and '0', we have assumed that they 
are always above the critical temperature Tc = 170 MeV 
even for large g. We note that the screening mass u i s 
nonzero in QCD vacuum but has a value of 180 MeV [59| . 
In this case, the binding energy of J in the vacuum is 
600~700 MeV. 
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FIG. 4: Charmonium dissociation temperatures in the QGP 
as functions of the QCD couphng constant g. 



Although the charmonium can be formed in the QGP 
at high temperature, it can still be dissociated by scat- 
tering with thermal partons. In the leading order (LO) 
pQCD, the charmonium breaks up by absorbing a ther- 
mal gluon, while in the next-to-lcading order (NLO) the 
dissociation is induced either by a quark or a gluon, and 
their invariant matrix elements arc given, respectively. 



by |61| 
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In the above, ki and k2 are, respectively, the momenta 
of incoming and outgoing thermal partons; tpip) is the 
wavefunction of charmonium with p = (ki — k2)/2; Nc 
is the number of colors; mc is the mass of thermal gluon 
and can be extracted from the lattice QCD [4l[; to$ is the 
mass of charmonium in the QGP; and nic = nic + (t/2/x 
is the mass of the constituent charm quark. With the 
screening mass fi = 0.18 GeV in the vacuum [HI], the 
latter has a value of mc = 1.85 GeV in the vacuum and 
is similar to the mass of D meson. The dissociation cross 
sections of charmonia are then obtained by integrating 
Eq. (fT9)) - ((2T|) over the phase space. 

The same pQCD formula can be used for charmo- 
nium dissociation by partons inside hadrons in the HG. 
It was found, however, that the charmonium is not heavy 
enough for pQCD to be applicable [12]. In the present 
study, we thus take the cross section for charmonium dis- 
sociation by a hadron to be proportional to its squared 
radius as in Ref . [l^l or given by that from a phenomeno- 
logical hadronic Lagrangian [63, H^l • We note that the 
effect of charmonium dissociation in the HG is negligible 
compared to that in the QGP due to the much smaller 
thermal decay width [l^, |43 . 

In terms of its dissociation cross section cr^'^^, the ther- 
mal decay width of a charmonium is given by 



r(r) 



E 



(2^ 



:V,,iik)n,{k,T)af^^k,T), (22) 



where i denotes the quarks and gluons in the QGP, and 
the baryons and mesons in the HG; rii is the number 
density of particle i in grand canonical ensemble; and 
Wrci is the relative velocity between charmonium and the 
particle. For the thermal width in the mixed phase, it 
is taken to be a linear combination of those in the QGP 
and the HG as following: 

r(T,) = / rQGP(T,) + (1 - /)rHG(T,), (23) 

where / is the fraction of QGP in the mixed phase. 
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FIG. 5: Thermal decay widths of charmonia in the QGP as a 
function of temperature for the QCD coupling constant g — 
1.87. 



The thermal decay widths of charmonia also depend 
both on the QCD coupling constant and the tempera- 
ture of QGP. In Fig. \E\ they are shown as functions of 
temperature for g = 1.87. It is seen that the thermal 
decay width of J/ diverges at the dissociation temper- 
ature T = 300 MeV, while those of Xc and of ip' become 
divergent at the critical temperature Tc = 170 MeV. An 
infinitely large thermal decay width implies that the par- 
ticles instantly reach their maximally allowed equilibrium 
value Np. Therefore, the J/ip abundance is not expected 
to reach this value at Tc, in contrast to that of the Xc and 
tjj' . We note that the value g = 1.87 is slightly larger than 
that used in our previous studies based on a schematic 
firecylinder model as a result of the viscous effect that is 
included in the present study. 



V. RESULTS 

Using the above described two-component model based 
on the schematic viscous hydrodynamics and taking into 
account the in-medium effects on charmonia, we can cal- 
culate the nuclear modification factor Raa of J/ip in 
heavy ion collisions according to 

Raa = — fxa ^ fip'{2S) — fb)Rpri + fbRb + Ricg, 

(24) 



where i?pri, Rb, and i?rcg are the nuclear modification 
factors for J/ip^s that are produced from primordial hard 
nucleon-nucleon scattering, the decay of bottom hadrons, 
and the regeneration in the QGP, respectively. In writ- 
ing the above expression, we have used the fact that all 
primordial Xc and tjj' are dissociated above the critical 
temperature Tc- For i?pri, it is calculated according to 



(25) 



where Tf is the freeze-out proper time. For Rh, it is 
taken to be one as a result of the expected conservation 
of total bottom and antibottom numbers. As to i?rcg, 
it is calculated from the ratio of the number of J ftp's 
obtained from solving Eq.(l5]) to the number of J ftp's from 
p-|-p collisions at same energy multiplied by the number 
of binary collisions in A-|-A collisions. 



Nuclear modification factor of Jftp at SPS and 
RHIC 



a: 



S 0.0 




400 



FIG. 6: Nuclear modification factor Raa of Jftp (solid line) 
as a function of the participant number A'part in Pb-|-Pb 
collisions at -y/sjvjv = 17.3 GeV at SPS (upper panel) and 
in Au4-Au collisions at y^SNN = 200 GeV at RHIC (lower 
panel). Dashed and dotted lines represent, respectively, 
the contributions to Jftp production from primordial hard 
nucleon-nucleon scattering and regeneration in the QGP. Ex- 
perimental data are from Refs. 

In Fig. ini we show the nuclear modification factor Raa 
of Jftp as a function of the participant number in Pb+Pb 
collisions at ^/sNN = 17.3 GeV at SPS and in Au+Au 
collisions at ^/snn — 200 GeV at RHIC. These results 
are obtained with the QCD coupling constant g = 1.87, 
which gives a good description of the experimental data 
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as shown by solid lines in the upper and lower panels. 
It is seen that the Raa of J /ip becomes smaller as the 
number of participants in the collision increases. Also 
shown in Fig. [S] are results from the primordial (dashed 
lines) and the regenerated J/i/' in the QGP (dotted lines), 
and they clearly indicate that the contribution from the 
primordial J /ip decreases and that from the regenerated 
ones increases as the collision energy increases. 

B. Nuclear modification factor of J/tp at LHC 



and shown by the dash-dotted line in Fig. [2]as a function 
of the participant number. This contribution is compa- 
rable to that from the regenerated J/ipi (dotted lines) in 
peripheral collisions and more important than the pri- 
mordial ones (dashed lines) in more central collisions. 
The upper and lower solid lines are the final Raa of J / ip 
obtained with the nuclear absorption cross section of 
and 2.8 mb, respectively. It is seen that the difference 
between the results obtained with and without the nu- 
clear absorption is mainly in collisions of small number of 
participants as the primordial J/ -^s are mostly dissolved 
in central and semi-central collisions. 
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FIG. 7: Nuclear modification factor Raa of J/t/i as a function 
of tlie participant number A''part witfiout (upper panel) and 
with (lower panel) tlie sliadowing effect in Pb-|-Pb collisions 
at Y^ijvjv = 2.76 TeV at LHC. Upper and lower solid lines 
are the Raa of J/ip for the nuclear absorption cross sections 
of cTabs =0 and 2.8 mb, respectively. Dashed, dotted, dot- 
dashed lines denote, respectively, the contributions to J/tp 
production from primordial hard nucleon-nucleon scattering, 
regeneration in the QGP, and decay of bottom hadrons. 

In Fig. [71 we show the Raa of J/ip in Pb+Pb colli- 
sions at ^SNN = 2.76 GeV at LHC with (lower panel) 
and without the shadowing effect (upper panel). It is 
seen that the shadowing effect suppresses the produc- 
tion of charm pairs and consequently the regeneration 
of J/ip- In obtaining these results, we have included the 
contribution to J/ ip production from the decay of bottom 
hadrons, which becomes non- negligible at LHC [sollioj. 
by assuming that the Raa of J/ip from the decay of bot- 
tom hadrons is independent of the centrality as indicated 
by the measured data from the CMS Collaboration [15 1 




5 10 15 

(GeV/c) 



20 



FIG. 8: Transverse momentum pr spectrums of J/tp from 
p + p annihilation at ^/sNN = 1-96 TeV by the CDF Collabo- 
ration at Fermi Lab (filled squares and solid line) [2^ and of 
regenerated J/^ (dashed line) in central Pb-|-Pb collisions at 
= 2.76 TeV at LHC. 



To compare the results from our model with the exper- 
imental data from LHC [lH , which have a transverse 
momentum cut px > 6.5 GeV for the measured J/tp, we 
note that the fraction of produced J/tjj's with transverse 
momentum larger than 6.5 GeV from the decay of bot- 
tom hadrons is about 21% in p -I- p annihilation at 1.96 
TeV from the CDF Collaboration at the Fermi Lab 
Parameterizing the latter by [1 -I- {px / 4:.lGeV)'^]~^'^ as 
shown by the solid line in Fig. [51 we obtain that the 
fraction of J/ip^s with transverse momentum larger than 
6.5 GeV is 3%. This is significantly larger than that 
from the regeneration contribution in Pb-|-Pb collisions, 
which is only 0.17%, as shown by the dashed line that is 
obtained from the two-component model but is arbitrar- 
ily normalized. It was first pointed out in Ref. [65j that 
limiting the J/ip transverse momentum to 6.5 GeV sup- 
presses the contribution from the regenerated J/tp- For 
J ftp's of high transverse momenta, their nuclear modifi- 
cation factor Raa can thus be calculated by multiplying 
the last term in Eq. (|M)) by the percentage of regener- 
ated J ftp's with transverse momentum larger than 6.5 
GeV divided by the percentage of primordial J ftp's with 
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the same range of transverse momenta, which is 0.12 in 
central Pb+Pb colhsions at ^/snn = 2.76 TeV. 



without shadowing 




a: 



^ 0.0 



with shadowing 




400 



FIG. 9: Nuclear modification factor Raa of J/tp with trans- 
verse momentum larger than 6.5 GeV versus the number of 
participants without (upper panel) and with (lower panel) the 
shadowing effect in Pb+Pb collisions at yijvjv = 2.76 TeV at 
LHC. Dashed, dotted and dot-dashed lines represent, respec- 
tively, the contributions to J/tp production from primordial 
hard nucleon-nucleon scattering, regeneration from QGP, and 
decay of bottom hadrons, and the solid line is the sum of them, 
without nuclear absorption. Upper and lower solid lines are 
the Raa obtained with the nuclear absorption cross section 
of and 28 mb, respectively. 

Since it takes time for an initially produced cc pair 
to form a charmonium, which depends on the charmo- 
nium radius and the relative velocity between charm and 
anticharm quark [6^, [g^ , thermal dissociation of char- 
monia is thus delayed until charmonia are formed. Since 
the J/ip formation time increases with its transverse mo- 
mentum as a result of time dilation, this effect becomes 
more important for J/ijj^s of high transverse momenta 
that are measured in experiments at LHC. In this study, 
we treat the formation time as a free parameter to fit the 
experimental data. Using the formation time of 0.5 fm/c, 
which corresponds to 1.4 fm/c in the firccylinder frame 
based on the average of the J/ijj transverse momenta that 
are above 6.5 GeV, our results for the the nuclear mod- 
ification factor Raa of J/ip with transverse momentum 
larger than 6.5 GeV as a function of the participant num- 
ber are shown in Fig. [S] without (upper panel) and with 
(lower panel) the shadowing effect as well as without (up- 
per solid curve) and with (lower solid curve) the nuclear 



absorption effect. It is seen that the results obtained 
without the shadowing and the nuclear absorption ef- 
fect describe reasonably the recent experimental results 
from the CMS Collaboration at LHC [l5| shown by solid 
squares. Also, it is interesting to see that the shoulder 
structure around iVpart = 100 in the measured Raa at 
LHC is roughly reproduced by our model. As suggested 
in our previous study on the J/ijj Raa in Au-|-Au col- 
lisions at RHIC, the sudden drop in its value at certain 
value of iVpart reflects the maximum temperature of the 
formed QGP that is above the dissociation temperature 
of J/ip, because the survival probability of J/ip is discon- 
tinuous at its dissociation temperature. Moreover, the 
fact that the shoulder seen at LHC occurs at a smaller 
number of participants than the value iVpart = 190 at 
RHIC is consistent with the expectation that the max- 
imum temperature of the QGP formed at LHC that is 
above the dissociation temperature of J/tp happens in 
more peripheral collisions than at RHIC. We note that 
the CMS Collaboration has also measured the fraction 
of J/ tp from the decay of bottom hadrons in Pb-|-Pb col- 
lisions at y^sJiN = 2.76 TeV and found that its nuclear 
modification factor is about 0.37, which corresponds to 
Rb in Eq. (p4|) . and is almost independent of the central- 
ity [H. 
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FIG. 10: Ratio Rep of the Raa of J/ip with transverse mo- 
mentum larger than 6.5 GeV in a given centrality to that in 
the peripheral collision versus the centrality of Pb-|-Pb colli- 
sions at y^SNN = 2.76 TeV at LHC. Experimental data shown 
by solid squares are from Ref. [l^ . 



Our results can further be compared with the experi- 
mental data from the ATLAS collaboration at LHC [l3| 
on the centrality dependence of the ratio i?cp of the Raa 
of J/iJj in a collision of certain centrality to that in the 
peripheral collision. For this purpose, we determine the 
centrality of a collision using the Glauber model as fol- 
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lows [68 



Centrality(6) jg''^ ' 



total incl 



AB 



/g°°27rfe'd6'|l - 



AB 



(26) 



where the numerator is the inelastic cross section of nu- 
clei A and B with the impact parameters between and 
b, and the denominator is the total inelastic cross section 
of the two nuclei; and ct;^^ is the inelastic cross section 
of a p+p collision at the same collision energy. In Fig. 
[TOl we show the calculated centrality dependence of Rep 
in Pb+Pb cohisions at = 2.76 TeV at LHC with 

the uncertainty of the reference point, i.e., the Raa of 
J/ip in the peripheral collision, shown as dashed lines. It 
is seen that results from our model calculations can re- 
produce the measured Rep of J/ip^s with high transverse 
momenta, and the shadowing and the nuclear absorption 
effect do not make significant difference in the i?cp of 



VI. SUMMARY 

Modeling the evolution of the hot dense matter pro- 
duced in relativistic heavy-ion collisions by a schematic 
viscous hydrodynamics, we have extended the two- 
component model, that was previously used to describe 
J/ip production in heavy- ion collisions at RHIC, to those 
at SPS and LHC. As in our previous studies, we have 
included the effect due to absorption by the cold nuclear 
matter on the primordially produced charmonia from ini- 
tial nucleon-nucleon hard scattering, the dissociation of 
survived charmonia in the produced hot dense matter, 
and the regeneration of chamronia from charm and an- 
ticharm quarks in the quark-gluon plasma. For heavy-ion 
collisions at LHC, we have further included the shadow- 
ing effect in the initial cold nuclei. We have also taken 
into account the medium effects on the properties of the 
charmonia and their dissociation cross sections by us- 
ing the screened Cornell potential model and the NLO 
pQCD. With the same quasiparticle model for the equa- 



tion of state of the QGP and the resonance gas model 
for that of the HG as used before, we have obtained a 
lower initial temperature than in our previous study to 
reach the same final entropy density as a result of the 
finite viscosity. Consequently, a slightly larger QCD cou- 
pling constant was needed to reproduce the measured 
centrality dependence of the nuclear modification fac- 
tor at RHIC. The calculated nuclear modification factor 
for heavy-ion collisions at the SPS was found to agree 
with the measured value as well. For both the SPS and 
RHIC, the contribution from the primordial charmonia 
was found to dominate, although the contribution from 
the regenerated ones increases from the SPS to RHIC. 
For heavy-ion collisions at LHC, the regenerated char- 
monia becomes most important in semi-central to cen- 
tral collisions as a result of the larger number of charm 
and anticharm quark pairs produced in higher energy col- 
lisions. Since the available experimental data from the 
LHC are for J/ip^s of transverse momentum px > 6.5 
GeV, we have further considered the contribution of J/ip 
production from the decay of bottom hadrons as its ef- 
fect increases with increasing J/'ijj transverse momentum 
and the effect due to the formation time of the J/ip- A 
reasonable agreement with the preliminary experimental 
data has been obtained if the shadowing and the nu- 
clear absorption effect is absent. However, a definitive 
conclusion can only be made after more refined exper- 
imental data becomes available. Furthermore, we have 
found the similar trend in the centrality dependence of 
the Raa of J/ ip at both RHIC and LHC that it decreases 
monotonously in peripheral collisions and then drops at 
a certain centrality as a result of the onset of an ini- 
tial QGP temperature higher than the J/tp dissociation 
temperature. Moreover, this takes place in less central 
collisions at LHC than at RHIC, indicating that the ini- 
tial temperature at the same centrality is higher at LHC 
than at RHIC. 
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